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ABSTRACT 



Context. Hydrodynamical cosmological simulations predict flows of the intergalactic medium along the radial vector of the voids, 

approximately in the direction of the infall of matter at the early stages of the galaxy formation. 

Aims. These flows might be detected by analysing the dependence of the warp amplitude on the inclination of the galaxies at the 

shells of the voids with respect to the radial vector of the voids. This analysis will be the topic of this paper. 

Methods. We develop a statistical method of analysing the correlation of the amplitude of the warp and the inclination of the galaxy 

at the void surface. This is applied to a sample of 97 edge-on galaxies from the Sloan Digital Sky Survey. Our results are compared 

with the theoretical expectations, which are also derived in this paper. 

Results. Our results allow us to reject the null hypothesis (i.e., the non-correlation of the warp amplitude and the inclination of the 

galaxy with respect to the void surface) at 94.4% C. L., which is not conclusive. The absence of the radial flows cannot be excluded at 

present, although we can put a constraint on the maximum average density of baryonic matter of the radial flows of (pj,) <~ 4f2/,p;.iit- 

Key words, intergalactic medium — galaxies: statistics — galaxies; kinematic and dynamics — large-scale structure of Universe 



1. Introduction 

Warps seem to be an almost universal structural feature in spiral 
galaxies. Indeed, most of the spiral galaxies for which we have 
relevant information on their structure (because they are edge 
on and nearby) present warps in their stellar and gas distribu- 
tions. Sanchez-Saavedra et al. (1990, 2003) and Reshetnikov & 
Combes (1998) show that nearly half of the spiral galaxies of 
selected samples are warped, and many of the rest might also 
be warped since warps in galaxies with low inclination are diffi- 
cult to detect. They are more clearly observed in the HI distribu- 
tion (see e.g. van der Kruit 2007 and references therein). Warps 
are also detected at about z = 1, even with a larger amplitude 
(Reshetnikov et al. 2002). 

Despite the compelling observational evidence of warps in 
the spiral discs, there is consensus on what could be the origin 
of this property of the galaxies. Nevertheless, it seems clear that 
warps should be produced by an interaction of the disc with an 
external element. In fact. Hunter & Toomre (1969) showed that 
in an isolated galaxy (without a dark matter halo), an initial warp 
would soon disappear and leave as its only trace a thickening of 
the edge of the disc. 

The number of ideas suggested to explain the origin of the 
warp in discs is vast. One explanation for the warps is gravita- 
tional tidal effects due to satellite galaxies. At least in the Milky 
Way galaxy, this explanation does not work with Magellanic 
Clouds as satellite (Hunter & Toomre 1969), and it is contro- 
versial whether it works in combination with the amplification 
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of the halo (as proposed by Weinberg 1998 and criticised by 
Garcia-Ruiz et al. 2002). Also, the intergalactic magnetic field 
has been suggested as the cause of galactic warps (Battaner et al. 
1990; Battaner et al. 1991; Battaner & Jimenez- Vicente 1998). 

Following the evidence that galaxies seem to be embedded 
in a massive dark matter halo, the interaction between the halo 
and the disc was explored. Ideas like 'dynamical friction' be- 
tween the disc and a spherical halo (Bertin & Mark 1980; Nelson 
& Tremaine 1995), a flattened halo misaligned with the disc 
(Toomre 1983; Dekel & Shlosman 1983; Sparke & Casertano 
1988; Kuijken 1991), or resonant interactions with a triaxial halo 
(Binney 1981) were explored. All these ideas, however, were re- 
jected when the dark matter halo was modelled correctly as a 
deformable mass of collisionless particles, rather than as a rigid 
body (Binney et al. 1998). Since a warp represents a misalign- 
ment of the disc's inner and outer angular momentum, Ostriker 
& Binney (1989) and Jiang & Binney (1999) proposed a model 
in which warps are generated through accretion of material into 
the halo with a misaligned spin that changes the major axis of 
the halo with respect to the disc and consequently produces a 
torque over the disc. There is a need for substantial accretion of 
low angular momentum material from the IGM into the galax- 
ies (Fraternali et al. 2007), and the direction of the net angular- 
momentum vector of the material that is currently being accreted 
should be constantly changing (Quinn & Binney 1992). 

Also based on infalling of material, but with a much weaker 
dependence on halo properties, some works (Mayor & Vigroux 
1981; Revaz & Pfenninger 2001; Lopez-Corredoira et al. 2002; 
Sanchez-Salcedo 2006) have proposed a mechanism for the for- 
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mation of the warp in terms of the infall of a very low den- 
sity intergalactic medium onto the disc without the dynamical 
intervention of an intermediate halo. Both S-type and U-type 
warps can be produced by this interaction (Lopez-Corredoira et 
al. 2002; Saha & Jog 2005). Even if there are other mechanisms 
able to produce warps, at least we know that the infall of material 
onto the disc will always produce warps. 

If the infall of material is relevant to the formation of the 
warp of the disc, the orientation of the galaxies within the cos- 
mological large-scale structure where they are embedded should 
have an effect on the formation of these features. There is grow- 
ing evidence that disc galaxies are not oriented randomly, but 
their angular momentum primarily point parallel to the filaments 
(or sheets) where they are located. In the supergalactic plane, 
there is a hint of an excess of galaxies whose angular momen- 
tum lie in this plane (Kashikawa & Okamura 1992; Navarro et 
al. 2004). Beyond the local universe, Trujillo et al. (2006) show 
at the 99.7% level that spiral galaxies located on the shells of 
the largest cosmic voids (r > 10 /i^'Mpc) have rotation axes that 
lie primarily on the void surface. Paz et al. (2008) point out that 
the angular momentum of flattened spheroidals in SDSS galax- 
ies tends to be perpendicular to the large-scale structure. These 
alignments are expected to be a consequence of the gain in an- 
gular momentum of the galaxies at the early stages of their for- 
mation, when both the baryonic component and the dark matter 
protohalo are suffering tidal torques from neighbouring fluctua- 
tions. Using N-body simulations, the alignments of the angular 
momentum of the haloes with the large-scale distribution have 
been also found (Porciani et al. 2002; Bailin & Steinmetz 2005; 
Brunino et al. 2007; Aragon-Calvo et al. 2007; Hahn et al. 2007; 
Paz et al. 2008). 

The aim of this paper is to check whether the orientation of 
the spiral galaxies in the void surfaces is related to the presence 
of a warp or not. In contrast to filaments (which are strongly 
affected by redshift-space distortion), large cosmological voids 
are a feature easy to characterise from the observational point 
of view. In addition, another important advantage of the void 
scheme is that (because of the radial growing of the voids) the 
vector joining the centre of the void with the galaxy position is a 
good approximation of the direction of the maximum compres- 
sion of the large-scale structure at that point. Consequently, the 
radial vector of the void at the galaxy position approximately 
represents the direction of the infall of matter at the early stages 
of the galaxy formation. At later epochs, however, most of the 
accretion of material in the galaxy is expected to be through the 
filaments (i.e. parallel to the void surface). According to Lopez- 
Corredoira et al. (2002), the infall of material should produce a 
correlation between the orientation of the galaxy and the ampli- 
tude and direction of the S -component, or the U-component or 
both of them. We want to check this hypothesis here. The aim 
of this paper is producing a method for analysing the relation- 
ship of the warp amplitude in galaxies with the inclination of 
the galaxy with respect to the line "centre of void"-galaxy (to 
check the early accretion of material). This method is then ap- 
plied to the edge-on galaxies and void catalogue used by Trujillo 
et al. (2006) for available images from Sloan Digital Sky Survey 
(SDSS) survey. The work presented here is an attempt to obser- 
vationally characterise the influence of the large-scale structure 
(and, consequently, the cosmic infall of material) on the forma- 
tion of the warps. Some works have previously dealt with no ran- 
dom orientations of warps on large scales (Battaner et al. 1991) 
or in the Local Group (Zurita & Battaner 1997), but not at the 
void shells. 



2. Definitions 

2.1. Warp amplitude 

To define a warp amplitude, we first rotate the galaxy to have 
the mean plane of the galaxy coincident with the constant dec- 
lination axis in the local plane of the sky (perpendicular to the 
line of sight). The position angle is calculated with an iterative 
method that fits the central part of galaxies (size of galaxy/2) to 
a straight line. This method uses the position angle from Trujillo 
et al. (2006) as starting point. The position angle was determined 
in this way to an accuracy of about 0.5 degrees. This error was 
adopted like that of the RMS using the mean least square method 
in the rotation procedure. We then have a right and a left part 
of the galaxy, each having its own warp. The right part is the 
one with a lower right ascension. To quantitatively estimate the 
warp amplitude we define the warp parameter W, on the left(l) 
or right(r) side of the galaxy, as 



W(r or V) - 



^""" xydx 



(1) 



where L(r or /) is the radius (left or right) of the disc within the 
limits in which the disc is visible cr = signal/noise > 3; cr 



- -Jcrfys + cr^jj, where osys is the systematic error and o-^td the 

standard deviation) and y is the height of the disc at position of 
pixel X,, being x > (Fig. [TJ. The y-values are obtained as the 
peaks of Gaussians fits in the light distribution perpendicular to 
the plane. As said, we only considered data with intensity greater 
than 3cr. An example of the result of our analysis is presented 
in Fig|4] where the warp curve is drawn only for those values 
with an error bar less than 0.5". This estimated error bar can 
be computed by scaling the standard deviation (Icr error) by the 
measured chi-squared value. Then, Wr will be positive for warp 
towards increasing declination, and vice versa for W/. A large 
warp can reach values of W = 5 x 10"^ and a barely perceptible 
warp W = 5 X lO"''. 

Expression ([TJ is adimensional, therefore the value of W only 
depends on the shape of the edge-on galaxy but not on the intrin- 
sic size or on the distance of a galaxy [neglecting the change of 
the factor (1 + z)'^ (i.e. cosmological dimming) in the surface 
brightness of the galaxies throughout our sample since most of 
them are at a similar z ~ 0.1]. For numerical purposes and work- 
ing in pixels, we use the discrete expression 



W(, 



'■ or I) 



(L, + L,y 



(2) 



where x,- > and Xf^ - L^i- or /)• The error in estimating W is 
dominated by the imperfect rotation step of the galaxy when it 
is rotated to make the major axis coincident with the x axis. As 
mentioned before, this rotation is performed with an error of 0.5 
degrees (i.e. about 0.008 radians). This introduces an error in W 
given by 



AW'(r or I) = 



J^^*' °'"x2 tan 0.008 rfx 



(Lr + L,)^ 



3.3x10" 



(3) 



Further details of this method of warp measurement are given in 
Guijarro et al. (2008). 

In a S-shape warped galaxy, W,- and Wi have the same sign. 
In a U-shaped galaxy, Wr and Wi have different signs. We define 
the variables S and U as 



S sWr + Wl, 



(4) 
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Fig. 1. Graphical representation of the warp measurement. 



V = Wr- Wl. 



(5) 



If the warp is of the type with integral-sign [S-warp, see Fig. 
|2lleft)], S will be different from zero, positive or negative, and 
U will be zero if it is perfectly symmetrical or has a low value 
if there is some asymmetry. Otherwise, if the warp is predomi- 
nantly cup-shaped [U-warp, see Fig.|2jright)], U will be different 
from zero and S zero or very low, since it expresses the degree 
of asymmetry with respect to a perfect U-shape. An L-warp will 
have|5| ^ |t/|. The combination ofS-warps and U-warps explain 
the asymmetry of the warps (Lopez-Corredoira et al. 2002; Saha 
& Jog 2006) and the values of 5 and U give us the degree of each 
component to the total warp. We assign a value of S and U to all 
the galaxies in our sample and make statistics with these num- 
bers, which quantify the S-component and the U-component. 

A serious difficulty arising in any observational study of 
warps is that companions, spiral arms, and other effects may 
mimic warps. The error introduced by a misidentification are 
difficult to evaluate. However, the images do not suggest that the 
warps are confused with spiral arms. On the other hand, the com- 
panions which are far away from the plane of the main galaxy 
are not confused with the warp, and if they were very close to the 
galactic outskirts, the galaxy would be removed from our list. 



2.2. Inclination of tine gaiaxy witii respect to tine centre of tine 
void 

For each galaxy, given its position angle and the position with re- 
spect to the centre of the void (see Trujillo et al. 2006 for details), 
we calculated the inclination of the rotation axis with respect to 
the line "centre of void"-galaxy. The sense of the rotation axis 
makes the "right" warp positive, that is, toward increasing dec- 



lination. And the inclination / is defined positive (between and 
n) if the line "centre of void"-galaxy is to the right (decreasing 
position angle) of the rotation axis or negative (between and 
-n) otherwise. Figure |2] illustrates this. 

The error in this inclination stems from the error on the dis- 
tance to the galaxies in Trujillo et al. (2006) sample. Due to the 
intrinsic motion of the galaxies away from the Hubble flow, this 
error is estimated to be around 4 /z 'Mpc, and the error in the 
distance to the centre of the void, around 2 /z^'Mpc. Taking into 
account that the average distance of the galaxies to the centre is 
» 12 /j"'Mpc. This leads to an average error of » 14°. Since 
these errors are statistical and not systematic, they will not affect 
the average signal that we find in the data, but will only decrease 
the signal-to-noise ratio. 

3. Theoretical predictions from a model of warp 
formation in terms of accretion of IGM onto the 
galactic disc 

3.1. Warp dependence on inciination of tine gaiaxy 

The predictions of the model with accretion of intergalactic 
medium (IGM) onto the disc for an average Milky Way-like 
galaxy is given in Lopez-Corredoira et al. (2002, Fig. 11). A 
fit of the curves in that figure gives the theoretical values 5, and 
U,: 



S-component: 



S,{e{deg.], y; < < 90) * 5o(v) 



(6) 



xsin(4.99- 1.266» + 0.16402 -0.003240^ -I- 1.99 X lO^V), 

SMdeg.], v; 90 < 6» < 180) = -5,(180 - 6, v); 

S,{e[deg.], v; -180 < < 0) = 5,(61 -H 180, v) 

U-component: 

U,{e[deg.],v\Q< 9<90)x t/()(v)cos(1.16 (7) 

-0.2950 -H 0.07966*2 - 0.0017361^ + 1.09 x 10" V), 

UMdeg.],v;90<9< 180) = -t/,(180 - 0, v); 

UMdeg.l y; -180 < < 0) = U,{-0, v), 

where 6 in these expressions is the direction of the IGM wind 
with respect to the rotation axis of the galaxy, and v the relative 
velocity of the wind. Together, Uo and 5o represent the maxi- 
mum amplitude of the S , and f/, that we calculate in the follow- 
ing sections. 

Assuming there is a wind flowing radially outwards in the 
void with velocity vT = 200 km/s (details will be given in 
Betancort-Rijo & Trujillo 2008), we must add the dispersion 
of velocities of the galaxies: o-j = 215 km/s, 0-2 - 209 km/s 
(Betancort-Rijo & Trujillo 2008) in the radial velocity vi (the 
projection of the velocity into the radial direction of the void) 
and the perpendicular component V2 with respect to the radial di- 
rection of the void with angular azimuth (p. To obtain these num- 
bers, Betancort-Rijo & Trujillo used the linear theory of grow- 
ing fluctuations in the large-scale structure. They computed the 
r.m.s. of the corresponding components of the velocity of mass 
particles on the surface of a void of lO/z 'Mpc with respect to its 
centre of mass. These numbers agree within a few per cent with 
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Fig. 2. Left: Graphical representation of a perfect S-warp (S + Q, U - 0). Right: Graphical representation of a perfect U-warp 
iU i^ Q, S - 0). North (higher declination) is up, south is down, "i" stands for the inclination between the Une "centre of void"- 
galaxy and the rotation axis of the galaxy. 



the numbers found in numerical simulations. Hence, the average 
warps are given by 



1 r*^ r*co r*co 

57(0 = - dcf> dv2 dviP(vi,V2)S,(e,v), 

^ Jo Jo J-oo 

lT,(i)^- f dcf, f dv2 f dviP(vuV2)Ur(0,v), 

^ Jo Jo J-oo 



6 - cos 



Pivuvi) 



vi cos ; + V2 sm ; cos < 



V 



V2 



V27ro-i 



■ exp 



(vi - ^f 
2cr? 



„2 \ 



2cr j 



= ^p(+^': 



(8) 

(9) 

(10) 

(11) 
(12) 



where Si{0,v) oc y^, Ut{0,v) oc v^ because the warp amplitude 
is proportional, both in the S-shape and U-shape, to v^ (Lopez- 
Corredoira et al. 2002; Eqs. (39), (45)). When we make all these 
calculations, we get a dependence of S ,{i) and U,(i) on / which 
is close to cos / (see Fig.O, although closer in the case of U,{i) 
than in the case of 5,(/). Therefore, from now onwards, we will 
consider as a first-order approximation [here we include the am- 
plitude resulting from the calculation with expressions ^, (|9]l 
approximately]: 



5r(0~ 0.985 COS!, 

TTrii) ~ l.lSC/oCOSi. 



(13) 
(14) 



To facilitate the comparison of our theory with the data, we 
estimate the correlations of 5/ and Ut [full expressions (O, ©] 
with the function cos /. 



(5, cos /) - (5,)<cos /) - 0.535 0, 
(L^cos i) - (U,){cos i) = O.eOUo. 



(15) 
(16) 



3.2. Amplitude of the S-component: So 

From Lopez-Corredoira et al. (2002, Figs. 10, 11, Eq. (39)), we 
can derive roughly that the maximum height y of the m = 1 com- 
ponent of warp of a Milky Way-like galaxy and baryonic mean 
density of the intergalactic medium pi, (roughly the average den- 
sity of the IGM flows radially ejected from the void to produce 
the observed effect): 

y = 2.8 X 10"^~^(km/s)pfc(kg/m^) exp [0.43 x(kpc)] kpc. (17) 

With the definition given in Eq. ([T]|, vi" - 200 km/s, multiplying 
by a factor l/n for averaging the integration of the line of nodes 
over all the angles, the maximum amplitude is 



|W(,o,/)|(0-O) = 



2.1 X 10>/,(kg/m") 



L(kpc)3 
X [(exp [0.43L(kpc)][L(kpc) - 2.33]) + 2.33] . 



(18) 



The size (semiaxis length) of a Milky Way-Uke galaxy is approx- 
imately L = 15 kpc. With this number. 



|W(,„,,)l(« = 0)~5xlO^Vi(kg/m^). 



That is, due to 5 = 2W, 



So ~ 10"pb{kg/m'). 



(19) 



(20) 



We must bear in mind that this is only an estimation of the order 
of magnitude, because not all the galaxies are like the Milky 
Way (although this is a reasonably good approximation for the 
galaxies in our sample). To give an example, with an average 
intergalactic medium density given by (p/,) = ^bPcrti ~ 8 x 10"^^ 
kg/m^ (taking O^ = 0.042; Spergel et al. 2007), we get a S- 
component with 5() ~ 8 x 10"''. 

3.3. Amplitude of the U-component: Uo 

Similarly, from Lopez-Corredoira et al. (2002, Figs. 10, 11, Eq. 
(45)), we derive roughly that the maximum height y of the m - 
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30 60 90 120 150 180 
ior(180+i)[deg.] 

Fig. 3. Dependence of S and U on / predicted by the theory. Solid line: expected trend if the warps were produced by intergalactic 
winds flowing radially outwards in the void, according to expressions ^ and ^, normalized to a maximum height of one. Dashed 
line: cos /. 



component of the warp of a Milky Way-like galaxy with IGM 
baryonic mean density p/, is 



<17tt-2 



y = 8.2 X 10"~ (km/s)p/,(kg/mOexp[0.38 Jc(kpc)] kpc. (21) 

With the definition given in ([B, vi" = 200 km/s, multiplying by a 
factor 2/n for averaging the integration of the line of nodes over 
the all angles, the maximum amplitude is 



\W^r or l)\(0 ^ 0) ^ 



7.0 X W'pbikg/m') 



L(kpc)3 

X [(exp [0.38L(kpc)][L(kpc) - 2.63]) + 2.63] . 
Wifli L = 15 kpc, 

|W(, or i)\(e = 0) ~ 0.8 X lO^'-phikg/m'); 
that is, due to U ^ 2\W\, 
Uo~L6xlO'-^Pb(kg/m'). 



(22) 



(23) 



(24) 



Again, we call that this is a rough estimation with Milky Way- 
like galaxies. With an average intergalactic medium density 
given by {pb) = ^bPcrit ~ 8 x 10"^** kg/m^ (taking Q.h = 0.042; 
Spergel et al. 2007), we get a U-component with L'^o = 1.3x10"^. 



4. Data and analysis 

We used the data of the SDSS-DR3 (3rd. data release) that have 
already been used in Trujillo et al. (2006). These data are their 
edge-on (inclination larger than 78°) galaxies, which are within 
the shells rn,id < r < rmid + 4 /z^'Mpc surrounding the largest 
voids, where r,,o,v/ > 10 /i^'Mpc is its radius. The lower the in- 
clination of the galaxy, the greater the thickness of the projected 
disc and, consequently, the greater the error in the determina- 
tion of the centroid of y{x). In the worst case (78°), the thick- 
ness of the projected disc is comparable to its intrinsic thickness 



(Dalcanton & Bernstein 2002), so the error in the warp ampli- 
tude is not significantly increased with respect to a 90° inclina- 
tion galaxy. 

The voids were located using maximal spheres empty of 
galaxies with magnitude over -19.31 - 5\ogh (Hq - lOOh 
km/s/Mpc), and they were found by means of the HB void finder 
(Patiri et al. 2006). From the SDSS available public data, we 
used the filter "r" images. In total we have 1 14 galaxies. 

For seventeen galaxies there were difficulties measuring the 
warp amplitude (for instance, due to the proximity of a star in the 
field or interaction with other galaxies), so there remain N -91 
galaxies with which we carried out the statistics (Table [TJ. In 
Fig.m we show three examples. 

Some images of warped galaxies could be the subject of al- 
ternative interpretations. For instance, considering the isophote 
maps, warp curve, and image in the central panel of Fig, H] a 
feature is found at ;c = -12, y = 5, either a companion galaxy 
or an inteloper, which could produce/modify the warp curve. 
However, we find that the warp at this galactocentric radius is 
real, as directly deduced from a detailed study of the isophote 
maps. 

If we plot their values of S and U vs. i, we get the results 
of Fig. ID There are sHght trends in S{i) [{S i(°)) - {S){i{°)) = 
(0.3° ± 6.1°) X 10-3] and in U{i) {{U i{°)) - {UXiC)} = (-8.6° ± 
5.6°) X 10^3]. The errors in the correlations are calculated as 
a-sCTil y/N and cruo-il Va^; where crs, cru and (Ti are the r.m.s. 
of the values of 5 , U, and /. 

The scattering of Fig.|5]may be for several reasons. For ex- 
ample, 

1 . The measurement of the warp amplitude has errors. 

2. Several mechanisms produce warps. The accretion produces 
signal and noise, while the other mechanisms only produce 
noise. 

3. Different masses of the galaxies not taken into account in our 
model, as said in ^3.2| and [3.3l l. 

4. The wind's other components apart from the radial one in- 
troduce scattering 
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Contour rrop ond S Worp Curve 



ontour mop ond No Worp Curve 




Fig, 4, Warp curves and contour maps of three selected galaxies (up); and 5-filters combined SDSS images of them (down), 50"x50". 
The isophotes are equidistant (in units of n*3cr equiv. to a step of +0.75 mag/arcsec^) starting at a level of about 3 cr above the 
sky background. Left: R.A. = 234.8898°, 5 = +3.2352° (J2000), S-warped. Centre: R.A. = 131.4264°, 6 = +51.2122° (J2000), 
U-warped downwards. Right: R.A. = 226.3795°, 6 = +57.4142° (J2000), negligible warp. 



The error of the inclination with respect to the void also in- 
troduces scattering. 



All these sources of contamination introduce an increase in 
the scattering but not a systematic error We can consider the 
error ofeach individual point in Fig. |5]as cTj - 1.7 x 10"^ for the 
S-component and ctu = 1.6 x 10"^ for the U-component (rm.s. 
in Fig. |5]). Our mission is to extract the statistical information 
hidden behind these clouds of points. 



4. 1 . Checking the null-hypothesis 

To check whether the null hypothesis is compatible with the ob- 
served U-components distribution, we computed the probability, 
P(/,), based on the binomial distribution of finding no more than 
Hq galaxies with U < Q and |/| < /,, and no more than n| galax- 
ies with U > and |/| > 180 - /», assuming that there is no 
correlation between / and U (i.e., the null-hypothesis). With this 
assumption, the probability that f/ > is 1/2 for any value of /, 
and using the binomial distribution for n|, «q, we find 



Pii.) 



1=0 ^ 'J\i=0 ^ 



( K 



m\ 



i-(mo+mi) 



(25) 



where mo is the number of galaxies with |/| < /»; nii is the number 
of galaxies with |/| > (180 - /,). To determine /», we assume that 
the signal is proportional to cos /, as determined previously in our 
model [see Eq. (fT4l i1. We assume that the galaxies are uniformly 



distributed in / (the small alignment reported by Trujillo et al. 
2006 is not very relevant to this purpose), then 



{U)(U) = ■ oc (1 + cos it). 

J^' sin i di 



(26) 



The r.m.s. is proportional to the number of galaxies within i < /,, 
which for a given sample and again assuming isotropy, has the 
proportionality 



(Tjjiit) oc (1 - COS /») 



-1/2 



(27) 



Thus, according to our model, the value of /, that maximizes 
the signal-to -noise ratio, ■^, is given for ;, = 60°. If instead 
of the approximate cosine dependence of Eq. ( ll4t . we took the 
exact calculation of Eq. (|9]l, the value that maximizes the sig- 
nal to noise ratio would be /, = 73° for U-component; for S- 
component, it would be i, = 75°. We assume as a good approxi- 
mation the cosine dependence, so we use /, = 60°. 

With this value of /, = 60°, the probability that our data 
are compatible with the null hypothesis is P = 0.056; that is, 
the null hypothesis is excluded within 94.4% C.L. If we took 
it - 45°, we would get P = 0.0043 (rejection of non-correlation 
within 99.57% C.L.), but this value of /, is not justified, a priori; 
therefore, the statistical significance must be less than this. For 
a higher value of /,, we also get rejection of the null hypothesis. 
For i, - 75° we get rejection within 95.7% C.L. 

If we do the same calculation for the S vs. / data, we 
find that the probability of null hypothesis cannot be rejected 
{P(i, - 60°) - 0.13). We also checked the null hypothesis 
with the Spearman rank correlation coefficient. This test gives 
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Fig. 5. Dependence of S and U on the inclination / in the observational data (stars). The squares with error bars represent the average 
in bins of / of 30 degrees. 



higher probabilities of a null correlation: P ■ 
P = 0.739 for 5(0. 



0.154 for t/(0 and 



4.2. Most likely values of the signal 

If we assume that Eq. ( fT4] i with some positive Uo applies, the 
mean value of U for / < /», U, is given by Eq. 



U = 0.88f/o. 



(28) 



The most likely value of U can be estimated from the data as 
follows. 

We assume that the distribution of the probabilities of a value 
of U, PiU), is Gaussian, centred at U with r.m.s. cru. For galax- 
ies with / < /,, the probability that f/ < is 



kJ — cx 



I.e., 



P{U)dU = ^ 



w; 



'^""li'-'H^. 



and the probability that t/ > is 



1 



I P(U)dU = - + w. 

Jo 2 



(30) 



(31) 



Identically, for / > 180-/, the probability that U <Ois ^+w, and 
the probability that U > is j - a>. The probability, F{lij), that 
the number of galaxies with / < /, and U < 0, k, is smaller than 
or equal to the observed value n^ (from a total of mo galaxies 
with / < it) and that, for ; > 180 - ;, the number of galaxies 
with U > Q, j, is lower than or equal to the observed value, «| 
(from a total of mi galaxies with i > 180 - /,), is given by the 
multiplication of the probabilities of both events: 



Fico)^ 



( "o 

Z 

U={) 



mo 
k 



(1/2 + «)""'"*( 1/2-w)* 



(32) 



Z 



m\ 
J 



il/2 + cj)""--'(l/2-ajy 



Thus, the mean value U is obtained from F[a)(U)] - 0.5, and the 
maximum and minimum values within 95% C.L. would be re- 
spectively T7+, TJ- derived from F[a>(U+)] = 0.95, F[w(I7_)] = 
0.05 respectively. Note that coiU) has been defined so that it must 
be positive if our model applies. A negative value of «([/) should 
be interpreted as evidence against it. 

_ With our data and_/, - 60°, the values are: U - 6.7 x 10""^, 
U+ = 14.5 X 10"^ U- = -0.3 X 10^4 Using our model, we 
can use these numbers to put a constraint on the IGM density. 



(29) From Eqs. ^ and (l24|, we find that (pt) = 4.8^^^ x 10 



-26 



kg/m"* (95% C.L.). This allows us to put an upper limit on the 
IGM density but not a minimum. The same calculation with 5- 
component gives a tighter constrain: 5 = 5.1x10""^, 5+ = 13.6x 
10-4, J_ = -2.3 X 10-4; <pi) = 6.7+iij^4 X 10-27 kg/m^ (95% 
C.L.). 



The correlation of S with a cosine function (approximately 
the expected shape theoretically) is {S cos /) - {S )(cos /) = 
(-0.004 ± 0.091) X 10-^ The correlation of U with U, (also ap- 
proximated to be a cosine function) is (U cos /) - (f/)(cos /) = 
(0.1 16 ±0.083) X lO""'. We can get a better constraint for the max- 
imum density from these correlations: with the S -component 
measurement and the expressions ( fTsT i and (|20] |. we find that 
(pt) <~ 3 X 10-27 kg/m^ * 4Q.bPcni (95% C.L.= 2cr); with the 
U-component measurement and the expressions ( fT6l ) and jTM . 
we find that <p/,) < 3 ~ lO^^e kg/m^ ^ 37Q.bPcri, (95% C.L.). 

Therefore, summarising the contents of this section, we re- 
ject the null hypothesis (i.e., the inclination of galaxies and the 
amplitude of the warp are not related to each other) at 94.4% 
C.L. Using our model, we can estimate the average density of 

the radial flow from the void to be O-M ^ookm/s ) ^bPcrit- 
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5. Conclusions 

Cosmological hydrodynamical simulations predict flows of IGM 
along the radial vector of the void. This radial direction is ap- 
proximately the same as the infall of matter in the early stages of 
the galaxy formation at the shells of the void. One way to search 
for the effect of this IGM flow in these shells is to measure the 
dependence of the warp amplitude on their galaxies as a func- 
tion of their inclination with respect to the radial vector of the 
void. In this paper, we have developed a method to measure that 
effect, and we made a first attempt to find this effect. The signal 
found in the U component of the warp (the null hypothesis is re- 
jected at 94.4% C.L.) gives some hint that such an effect might 
exist. This result is not conclusive (5.6% is not a very negligi- 
ble probability) and the absence of the radial flows cannot be 
excluded at present. If the IGM radial flows in the radial direc- 
tion of the voids exist, their baryonic matter density should be 
ipb) <~ 3 X 10"^^ kg/m"' - 4-Q.hPcrit- This density would increase 
inversely proportional to the square of the mean flow velocity if 
its value differs from 200 km/s. There is also the possibility that 
the accretion of material have diff'erent initial velocities than the 
radial direction of the void. 

There may be other mechanisms of warp formation different 
to the accretion onto the disc, but they would produce noise in 
the correlation if they have nothing to do with the IGM accre- 
tion. If the correlation of S-component amplitude and inclina- 
tion were observed, although it would be an argument in favour 
of Lopez-Corredoira et al. (2002) theory, it would not be to- 
tally conclusive because there might be alternative explanations 
for the correlation. The mechanism of accretion into the halo 
(Ostriker & Binney 1989; Jiang & Binney 1999) rather than onto 
the disc might possibly explain the correlation. There might be a 
relationship between warps and filaments associated to the void 
produced by primordial magnetic fields, or the frozen magnetic 
fields were aligned with the filaments (Florido & Battaner 1997), 
if the magnetic fields are also responsible for the warp formation 
(Battaner et al. 1990; Battaner et al. 1991; Battaner & Jimenez- 
Vicente 1998). However, these theories do not explain the U- 
component (the asymmetry of the S-warps), which are clearly 
observed in many galaxies (e.g., Reshetnikov & Combes 1998; 
Sanchez-Saavedra et al. 2003). The trend in the correlation of 
the U-component with the inclination of the galaxy obtained in 
this paper, if confirmed with higher statistical significance, could 
be taken as confirmation that the mechanism of IGM accretion 
onto the disc produces warps. The application of the method pre- 
sented in this paper to galaxy samples with more objects and/or 
better measurements of the warp amplitude is expected to give 
more accurate results. 

Acknowledgments: Thanks are given to the anonymous ref- 
eree for helpful comments, and to Joly Adams (language ed- 
itor of A&A) for proof-reading this paper Funding for the 
creation and distribution of the SDSS Archive has been pro- 
vided by the Alfred P. Sloan Foundation, the Participating 
Institutions, the National Aeronautics and Space Administration, 
the National Science Foundation, the U.S. Department of 
Energy, the Japanese Monbukagakusho, and the Max Planck 
Society. The SDSS Web site is http://www.sdss.org/ The SDSS 
is managed by the Astrophysical Research Consortium (ARC) 
for the Participating Institutions. The Participating Institutions 
are The University of Chicago, Fermilab, the Institute for 
Advanced Study, the Japan Participation Group, The Johns 
Hopkins University, the Korean Scientist Group, Los Alamos 
National Laboratory, the Max-Planck-Institute for Astronomy 



(MPIA), the Max-Planck-Institute for Astrophysics (MPA), New 
Mexico State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, the United States Naval 
Observatory, and the University of Washington. MLC was sup- 
ported by the Ramon y Cajal Programme of the Spanish Science 
Ministery. We thank the Spanish Science Ministery for support 
under grant AYA2007-67625-CO2-01. 



References 

Aragon-Calvo. M. A., van de Weygaert, R., Jones, B. J. T., & van der Hulst, J. 

M. 2007, ApJ, 655, L5 
Bailin, J., & Steinmetz, M. 2005, ApJ, 627, 647 
Battaner, E., Florido, E., Sanchez-Saavedra, M. L. 1990, A&A, 236, 1 
Battaner, E., Garrido, J. L., Sanchez-Saavedra, M. L., & Florido, E. 1991, A&A, 

251,402 
Battaner, E., & Jimenez- Vicente, J. 1998, A&A, 332, 809 
Betancort-Rijo, J., & Trujillo, I. 2008, in preparation 
Bertin, G., & Mark, J. W.-K. 1980, A&A, 88, 289 
Binney, J., 1981, MNRAS, 196, 455 

Binney, J., Jiang, I.-G., & Dutta, S. 1998, MNRAS, 297, 1237 
Brunino, R., Trujillo, I., Pearce, F. R., & Thomas, R A. 2007, MNRAS, 375, 184 
Dalcanton, J. J., & Bernstein, R. A. 2002, AJ, 124, 1328 
Dekel, A., & Shlosman, I. 1983, in: Intern, kinemat. & Dynam. galaxies (lAU 

Sym. 100), p. 187 
Florido, E., & Battaner, E. 1997, A&A, 327, 1 
Fraternali, P., Binney, J., Oosterloo, T, & Sancisi, R. 2007, New Astron. Rev. 

51,95 
Garci'a-Ruiz, I., Kuijken, K., & Dubinski, J. 2002, MNRAS, 337, 459 
Guijarro, A., Peletier, R., Battaner E., Jimenez- Vicente, J., de Grijs, R., & 

Florido, E. 2008, in preparation 
Hahn, O., Porciani, C, Carollo, C. M., & Dekel, A. 2007, MNRAS, 375, 489 
Hunter C, & Toomre, A. 1969, ApJ, 155, 747 
Jiang, I.-G., & Binney, J. 1999, MNRAS, 303, 7 
Kashikawa, N., & Okamura, S. 1992, PASJ, 44, 493 
Kuijken, K. 1991, ApJ, 376, 467 
Lopez-Corredoira, M., Betancort-Rijo, J. E., & Beckman, J. E. 2002, A&A 386, 

169 
Mayor, M., & Vigroux, L. 1981, A&A 98, 1 

Navarro, J. F, Abadi, M. G., & Steinmetz, M. 2004, ApJ, 613, L41 
Nelson, R. W., & Tremaine, S. 1995, MNRAS, 275, 897 
Ostriker E. C, & Binney, J. J. 1989, MNRAS, 237, 785 
Patiri, S., Cuesta, A. J., Prada, F, Betancort-Rijo, J., & Klypin, A. 2006, ApJ, 

652, L75 
Paz, D., Stasyszyn, F, & Padilla, N. 2008. larXiv:0804.4477 1 
Porciani, C, Dekel, A., & Hoffman, Y. 2002, MNRAS, 332, 339 
Quinn, T, & Binney, J. 1992, MNRAS, 255, 729 
Reshetnikov, V, Battaner E., Coinbes, F., & Jimenez- Vicente, J. 2002, A&A, 

382,513 
Reshetnikov V., & Coinbes F 1998, A&A, 337, 9 
Revaz Y., & Pfenninger D., 2001, in: Gas and Galaxy Evolution (ASP Conf. 

Proc. 240), J. E. Hibbard, M. Rupen, J. H. van Gorkom (eds.), San 

Francisco: Astronomical Society of the Pacific, p. 278 
Saha, K., & Jog, C. J., 2006, A&A, 446, 897 

Sanchez-Saavedra, M. L., Battaner, E, & Florido, E. 1990, MNRAS 246, 458 
Sanchez-Saavedra, M. L., Battaner, E., Guijarro, A., Lopez-Corredoira, M., & 

Castro-Rodri'guez, N. 2003, A&A 399, 457 
Sanchez-Salcedo, F 2006, MNRAS 365, 555 
Sparke, L. S., & Casertano, S. 1988, MNRAS, 234, 873 
Spergel, D. N., Bean, R., Dore, O., et al. 2007, ApJS 170, 377 
Toomre, A. 1983, in: Internal kinematics and dynamics of galaxies, D. Reidel 

Publishing Co., Dordrecht, p. 177 
Trujillo, I., Carretero, C, & Patiri, S. 2006, ApJ 640, LI 1 1 
van der Kruit, P C. 2007, A&A 466, 883 
Weinberg, M. D. 1998, MNRAS, 299, 499 
Zurita, A., & Battaner, E. 1997, A&A 322, 86 



Lopez-Corredoira et al.: Warp / voids 

Table 1. Amplitude of the warp (with the corresponding error), *-^°"'- tablelU 

in units of 10"^, and inclination with respect to the radial direc- r.a. Decl. W, Wi Err(Wr/i) iC) 

tion of the void of the used SDSS galaxies in this paper. 195.1273 0.4801 UA 03 03 -92.6 

195.7355 3.7370 -1.3 0.2 1.4 -146.4 

195.7591 -2.4564 -0.6 -0.3 1.1 105.5 

R.A.(°) Decl.(°) W,. Wi Err{Wr/i) iC) 197.1058 61.4598 0.6 -0.7 0.9 -126.0 

-63.1 197.8032 49.4153 0.0 -1.4 0.9 -99.3 

91.6 198.0996 -2.8867 1.6 -0.2 1.4 -80.0 

-55.6 198.5669 -2.0749 0.1 0.0 1.7 -60.9 

-51.1 198.9388 59.8330 0.8 0.5 1.1 -13.3 

-140.5 199.0415 65.9117 -4.3 0.2 2.3 -94.9 

-133.5 204.4937 5.7839 0.6 1.3 2.3 87.1 

-113.5 207.0692 46.7701 -1.2 -0.5 1.1 155.8 

79.1 207.7140 5.2116 -0.9 -0.4 1.3 128.0 

107.6 212.5481 58.7175 1.8 -0.5 1.9 -41.1 
-120.4 214.8475 2.8468 2.5 -0.7 1.5 139.9 

110.3 219.3559 0.4227 0.2 0.1 0.5 -100.3 
-77.7 220.0357 3.0825 -0.3 -0.1 0.3 68.2 
94.0 220.3983 62.7731 3.8 -0.4 2.8 -79.1 
-79.5 220.7501 61.6237 1.0 -0.5 1.6 119.9 
87.6 223.3199 49.1950 -0.2 1.5 1.5 163.5 

112.7 225.0443 43.1852 0.6 -0.3 1.6 148.9 
-72.8 226.3795 57.4142 1.3 0.3 1.2 99.4 
-113.0 226.8317 0.5954 -1.3 -1.1 1.1 -105.2 
99.0 227.2538 41.3454 -1.4 -0.5 1.6 71.3 
95.0 230.8046 39.8048 0.7 -1.2 1.5 86.2 
-144.9 233.1236 50.0678 -0.8 1.9 1.7 74.1 
-65.6 233.7380 58.4997 -0.3 0.0 0.8 48.6 
91.5 234.8898 3.2352 -2.9 -1.8 2.4 -117.9 
-80.5 236.9465 55.5457 0.6 -0.8 1.5 -28.5 
-77.2 238.5349 45.0029 0.0 0.3 1.3 -21.1 

36.5 239.2935 45.2622 0.3 0.0 0.7 88.7 
-129.4 241.3029 42.8759 -0.9 1.1 1.5 -61.0 

124.4 242.0378 50.1382 -0.1 -0.1 1.2 -76.3 
-156.8 243.6710 38.0918 -0.2 -0.9 1.6 -31.5 
4.0 244.8777 42.4485 1.3 -1.0 1.7 -20.4 

64.6 249.2328 44.0929 1.4 0.3 2.1 -36.1 
43.6 251.6845 42.7894 0.0 1.2 1.9 -56.0 
-147.2 254.6449 32.6595 0.0 0.6 0.7 -66.5 

100.3 257.1461 59.3407 -0.2 -0.6 0.7 -94.4 

100.4 258.1745 30.4143 -2.1 -0.2 3.2 112.7 
-101.6 324.0977 -6.4763 -1.9 0.0 0.8 108.4 
-128.1 341.0004 -0.9310 1.3 -0.4 3.0 -141.4 
107.2 354.9044 14.5762 -1.1 0.5 1.3 -76.9 
94.7 
32.3 
73.5 
-36.8 
80.4 
152.2 
-75.4 
56.5 
-91.0 
-113.9 
-109.8 
-117.5 
-34.1 
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65.4 
88.0 
81.2 
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